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Abstract 

We employ Coleman- Weinberg and Higgs potentials to implement inflation in non- 
super symmetric grand unified theories (GUTs) such as SU(5) and SO(10). To realize a 
scalar spectral index close to 0.96, as indicated by the most recent WMAP 5-year anlaysis, 
the energy scale of observable inflation turns out to be of order 10 16 GeV. This implies 
a GUT symmetry breaking scale of similar magnitude, and proton lifetime of order 10 34 - 

10 38 

years. In some SO(10) models with axion dark matter, the scalar leptoquark boson 
exchange leads to proton decay with a lifetime of order 10 -10 35 years. 

Introduction 

An inflationary scenario [lj may be termed successful if the following conditions are satisfied: 

1) The number of e-foldings is sufficient to resolve the horizon and flatness problems. 

2) The predicted temperature anisotropies, scalar spectral index n s , tensor-to-scalar ratio r, 
etc., are in agreement with the observations. 

3) Following inflation, there exists a mechanism for explaining the observed baryon asymmetry. 

4) The number density of superheavy magnetic monopoles in GUTs is suitably suppressed to 
agree with the upper bounds on their primordial abundance. 

5) The model should offer a plausible cold dark matter candidate. 

It was shown a long time ago [2l[3jH] that non-supersymmetric GUT inflation can be realized 
by employing a Coleman- Weinberg (CW) type potential for a scalar inflaton field which must be 
a GUT singlet. The singlet condition ensures that radiative corrections from the gauge sector 
do not spoil the desired inflationary potential. The monopole problem is resolved by requiring 



that the GUT symmetry breaking, associated with monopole production, takes place during 
inflation. A suitable coupling between the GUT symmetry breaking field £ and the inflaton 
field is thus needed. This coupling also plays a crucial role in generating an appropriate CW 
inflationary potential for <p. 

The most recent WMAP 5-year [5] analyses seem to favor a scalar spectral index centered 
around 0.96. Within the CW framework this favors an energy scale for inflation which is close 
to (1-2) xlO 16 GeV [H]. This result, as we will show here, has far reaching implications for 
the symmetry breaking scale of the underlying GUT. In our CW examples below, primarily 
based on SO(10), the superheavy gauge bosons which mediate proton decay are typically a 
factor 2-4 larger in mass than the energy scale of inflation preferred by WMAP5. To phrase 
things somewhat differently, while the most recent proton lifetime limits from Super-K (r(p — > 
vr°e + ) > 8.2 x 10 33 years [7j) require that M x > 4 x 10 15 GeV [8], the limits on M x from 
WMAP seem to favor a value which is up to a factor 10 larger. We thus estimate a proton 
lifetime of around 10 34 -10 38 years. Proton lifetime estimates of order 10 34 years are achieved 
in some SO(10) axion models which contain scalar leptoquarks (3,1,-2/3) with an intermediate 
mass of order 10 12 GeV. 

As an alternative to the CW potential, we will also discuss inflation with a Higgs potential 
[HI [I], as well as implications for proton decay. As we will see, such a model predicts values for 
the scalar spectral index n s , the tensor-to-scalar ratio r, the scale of vacuum energy and the 
proton lifetime that are quite similar to those predicted by the CW model. 

GUT Inflation with a Coleman- Weinberg (CW) Potential 

To simplify the presentation we begin with a brief summary based on SU(5), following references 
[2] and [10]. The model contains an SU(5) singlet (real) scalar <fi which develops a CW potential 
from its weak couplings to the adjoint and fundamental Higgs fields £ and H 5 . The tree-level 
scalar potential is given by 



The coefficients a, b, a and (3 are taken to be of order g 2 , so that most radiative corrections in 
the (E, H 5 ) sector can be neglected. We assume a somewhat smaller value for the coefficient 7 
and take < A« <C g 2 and Ai < max(A|, A|). 

Radiative corrections due to the couplings 4> 2 Tr £ 2 and 2 H\ H 5 induce a Coleman- Weinberg 
potential for 0, which is given by 



V(<f>,X,H 5 ) 



= 1 a ( Tr S 2 ) 2 + l - b Tr E 4 - a {HIE,) Tr E 2 + | (fljff 5 ) 2 
+ 7^S 2 iy 5 + ^ 4 -y0 2 TrE 2 + ^ 2 iytF 5 . 
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and A c = a + — b. The SU(5) symmetry is broken to SU(3) C x SXJ(2) L x U(l)y when E 

1 

acquires a vacuum expectation value (VEV) 



(£) = jia-diagl 1,1,1,-=— ) . Hi 



where 



2 2 ^2 i2 /r\ 

o- = — • (5) 
With E given in Eq. PJ, the (0, cr) sector of the effective potential can be written as 

V ■=£« 4 -£*V + ^ '(»-(^)+c)+V„. (6) 

Using Eq. (jSJ) and an appropriate choice of the normalization constant C, the effective potential 
for the inflaton field can be expressed in the standard form [10] 



V(0)=A0 4 



In 
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where M = (0) denotes the VEV of at the minimum. Note that the vacuum energy density 
at the origin is given by Vo = -^¥- such that V(0 = M) = 0, and the corresponding minimum 

m a is located at a = J^M. The mass of the superheavy gauge bosons X which mediate 
proton decay is given by 



5 g ct / 5A 2 ff 2 1/4 



Thus is estimated to be a factor 2-4 larger than the scale of vacuum energy during inflation. 
This is to be contrasted with the simplest supersymmetric hybrid inflation models in which the 
corresponding Mx can easily be 1-2 orders of magnitude larger than the scale of vacuum energy 
during inflation [TTj . 

The inflationary potential in Eq. (J7|) is typical for the new inflation scenario |12j . where 
inflation takes place near the maximum. However, as shown in reference [6], depending on the 
value of Vo, the inflaton can have small or large values compared to the Planck scale during 
observable inflation. In the latter case, observable inflation takes place near the minimum of 
V(0), and the model mimics chaotic inflation [T5j. Indeed, we will see from Table [TJ that in 
order to obtain a scalar spectral index close to 0.96, the energy scale for observable inflation is 
typically on the order of 10 16 GeV. 

So far, we have been discussing the case in which <fi < M during inflation. Alternatively, 
the inflaton may roll toward its minimum starting from values larger than the VEV, similar to 
simple models of chaotic inflation [U]. This is true in both of the inflationary models (CW and 
Higgs) that we consider here. For shorthand, we henceforth denote these regimes as the BV 
(below VEV) and AV (above VEV) branches. 
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Results for the CW Model 



The inflationary slow-roll parameters are defined as [2] 

1 {V'\ 2 fV"\ _ 2 (V'V"'\ 

(Here and below we use units mp = 1, where mp ~ 2.4 x 10 18 GeV is the reduced Planck 
mass.) The slow- roll approximation is valid as long as the conditions e <C 1 and rj <C 1 hold. 
In this case the spectral index n s , the tensor-to-scalar ratio r and the running of the spectral 
index 4^ are given by 



n s ~ 1 - 6e + 2r] (10) 
r ~ 16e (11) 

I6 eri _ 24e 2 - 2£ 2 . (12) 



(ilnfc 



The number of e- foldings after the comoving scale Iq = 27r//c has crossed the horizon is given 

where 0o is the value of the field when the scale corresponding to k exits the horizon, and 
<p e is the value of the field at the end of inflation. The value of <ft e is given by the condition 
2(H'(<f) e )/ H(<p e )) 2 = 1, which can be calculated from the Hamilton- Jacobi equation 



[H'm 2 - ^ 2 (0) = ~\V{4>) . (14) 
The amplitude of the curvature perturbation is given by 

1 v 3 / 2 

^^T,W\- (15) 

(Note that, for added precision, we include in our calculations the first order corrections in the 
slow-roll expansion for the quantities n s , r, and |16j.) 

To calculate the magnitude of A and the inflationary parameters, we use these standard 
equations above. The WMAP5 value for A^ is 4.91 x 1(T 5 for k = 0.002 Mpc" 1 [5j. The 
observable number of e-foldings corresponding to the same scale is 
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vl/4 



10 15 GeV 



+ - In 

3 



T 

± r 



10 9 GeV 



(16) 



For our calculations here and in the next section, we use the formula given in Eq. ( 1281) to 
obtain values for the reheat temperature T r . This formula is based on an SO(10) GUT model, 
and yields values on the order of 10 6 -10 7 GeV. In contrast, SU(5) considerations give T r ~ 10 9 
GeV. 
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y 1/4 (GeV) 


V^ ) 1/4 (GeV) 


A (10~ 14 ) 


M 


00 
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r 


T r (GeV) 


ttt(-10~ 3 ) 


2. x 10 15 


2.00 x 10 15 


3.7 


2.65 


0.11 


2.05 


0.9369 


0.000013 
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18.1 


0.9578 
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Table 1: Predicted values of various inflationary parameters using the Coleman- Weinberg po- 
tential given in Eq. ([7]). Here we show only those values which fall inside the WMAP5 2a 
bounds (see Fig. [TJ), and also omit values with very low reheat temperature T r (< 10 6 GeV). 
Note that unless otherwise specified, we use units of mp = 1. Although M, O and e carry 
transplanckian values, the vacuum energy scale during observable inflation is well below nip. 



Our predictions for the values of various parameters are displayed in Table [TJ As we see 
in Fig. [TJ both the BV and AV branches fall within the WMAP5 2a bounds for some range, 
and part of the BV branch lies inside the la bound. As discussed earlier, larger values of V^ 4 
result in observable inflation occurring closer to the VEV M. For the largest values of V^ A , the 
potential is effectively given by V = |m^(A0) 2 , where A0 = M — denotes the deviation of 

the field from the minimum and = 2\/~AM ( — J is the inflaton mass in the CW model. 

This well-known monomial model [13] predicts ~ 2 x 10 13 GeV, A0o — 2a/^o, n s — 1 — 
and r ~ 4(1 — n s ), corresponding to V(4>o) — (2 x 10 16 GeV) 4 . This is the region in which the 
two branches meet, i.e. both the BV and AV branches converge to quadratic inflation in the 
high-Vo limit. 

In the opposite limit with lower values of Vq, the BV branch reduces to new inflation models 
with an effective potential V = V (l — |A^|(-g) 4 ), where A^ = In (-g) , leading to n s ~ 1 — 

and r ~ y(l — n s )|A^ | (fj) 4 ~ 0. On the other hand, the AV branch is asymptotic to quartic 
inflation with an effective potential V = (^) 4 , resulting in n s ~ 1 — and r ~ y (1 — n s ). 

As previously indicated, Table [TJ shows that the energy scale of observable inflation in this 
model is ~ 10 16 GeV. It is instructive to see how n s changes with this quantity, and this is 
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0.92 0.94 0.96 0.98 1.00 1.02 1.04 

n s 

Figure 1: r vs. n s for the Higgs and Coleman- Weinberg models, shown together with the 
WMAP5 contours (68% and 95% confidence levels) [5] . In each model, inflation is allowed both 
below the VEV (BV) and above the VEV (AV). The BV and AV branches for the Higgs (CW) 
potential are shown in red (brown) and blue (green), respectively. The black circle corresponds 
to a quadratic potential with N = 52, n s = 0.961 and r = 0.152. 



displayed in Fig. [2j 



GUT Inflation with a Higgs Potential 

In this section we implement inflation by employing a Higgs potential given by [HI El I 

V{<P) = V -ii 2 <p 2 + \<p\ (17) 

It is useful to rewrite this in terms of the vacuum potential Vq and the vacuum expectation 
value M of the inflaton [HI I 



V (0) = v„ 



1-1*. 

M 



-I 2 



As mentioned earlier, and previously discussed in reference (9j, inflation may occur above 
or below the VEV M. We proceed in the same way as before, noting that Vq is not related to 
other parameters in a simple way as it was in the CW case. In this model, the inflaton mass is 
given by 



2 v / 2Vb 



M 



The results of this calculation are given in Table [2j For small values of Vo, the AV branch 
approaches quartic inflation with an effective potential V = (jg?) 4 , which is the same as in 
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Figure 2: n s vs. log lo [V A (0 o ) 1 ^ 4 /GeV] for the CW and Higgs models. The BV and AV branches 
for the Higgs (CW) potential are shown in red (brown) and blue (green), respectively. The 
black circle corresponds to a quadratic potential with N = 52 and n s = 0.961. 



the CW case apart from the slowly varying function A^. On the other hand, the BV branch 



by WMAP5. This leads to n s ~ 1 - ^ In (jj^J and r ~ 8(1 - n^e - ^ 1- "*). In contrast 

to the AV regime, the CW and Higgs models produce distinct predictions in the BV regime, 
which is apparent from Fig. [TJ For large Vo, both branches converge to quadratic inflation and 
to the CW solution. As in the CW model, observable inflation in this region takes place near 
the VEV M. Hence, it is not surprising that the two models yield very similar results. 

As shown in Fig. [TJ the Higgs model yields results compatible with WMAP5 for a wide range 
of values. Analogous to the CW case, BV inflation appears to be somewhat more favorable 
than AV inflation, although both branches fall within 2a of the central value for some range. 
In this model, the energy scale of observable inflation is again ~ (1-2) x 10 16 GeV, except for 
regions far outside the WMAP bounds. This can be seen in Table [2] for the most favored range 
of (r,n s ) values, as well as in Fig. [2j 

Note that radiative corrections to the inflationary potential, especially from Yukawa inter- 
actions, can modify the tree-level inflationary predictions in Table [2J This has recently been 
discussed for (p 2 and 4 potentials in reference [T9~] . 



approaches new inflation with an effective potential V — Vq ( 1 — 2 (-^j) 



in a region disfavored 
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Table 2: Predicted values of various inflationary parameters using the Higgs potential given in 
Eq. f|T8|) . analogous to Table [TJ 

Realistic SU(5) and SO(IO) Models 

It is well known that satisfactory gauge coupling unification is not achieved in the minimal 
non-supersymmetric SU(5) model |20j. However, this situation can be improved by introducing 
additional matter fields at a low energy scale. For example, following reference [21], if we 
introduce vectorlike fermions Q(3, 2, 1/6) + Q(3, 2, -1/6) and D c (3, 1, 1/3) + D c (3, 1, -1/3) 
with masses ~ 200 GeV, the gauge couplings unify at a scale Mx ~ 3.8 x 10 16 GeV, as shown 
in Fig. [3] Now using V^ 4 ~ Mx/2 and referring to Tabled], we obtain a corresponding value 
of the spectral index n s ~ 0.962 which is well inside the la bound of WMAP5 (see Fig. [[]). 

As mentioned earlier an inflationary scenario should provide, among other things, both a 
suitable cold dark matter candidate and an explanation of the observed baryon asymmetry. 
We will now consider a class of SO(10) GUT models in which these conditions are readily met 
[22J. In these models axions comprise the dark matter of the universe, and the observed baryon 
asymmetry arises via leptogenesis [23] [21] . Clearly, axion dark matter can also be introduced 
in the SU(5) model. 

The breaking of SO (10) to the Standard Model (SM) proceeds via an intermediate step 
SU(3) C x SU(2) L x SU(2) R x U(1) B _ L [25]. The U(1) PQ (Peccei-Quinn) symmetry and 
SU(2)^xU(1)b-l are both spontaneously broken at some scale Mb-l- Axion dark matter 
physics requires that Mb-l is around 10 11 -10 12 GeV [26]. Furthermore, from our earlier dis- 
cussion, satisfactory CW inflation requires that the SO(10) symmetry breaking scale be close 
to, indeed somewhat larger than 10 16 GeV. It is intriguing that the scale of B — L and axion 
symmetry breaking and the inflaton mass are of the same order of magnitude. This enables 
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Figure 3: Gauge coupling unification in the SU(5) model with additional fermions Q + Q + 
D c + D c at mass scale ~ 200 GeV. 



the inflaton to produce right-handed neutrinos whose subsequent decay produces the observed 
baryon asymmetry via non-thermal leptogenesis [24|. 

The fermion content of the SO(10) model consists of three SM families in the 16-dimensional 
spinor representations, as well as two fermion matter multiplets in the 10-plet representations 
[22J. These 10-plet fields are included in order to resolve the well-known axion domain wall 
problem [27], by ensuring that a residual, discrete PQ symmetry coincides with the center, Z4, 
of SO(10) [28]. Under U(1)pq, the fermion fields transform as follows: 

irf? — e W $l U = 1, 2, 3), 4o ] — e^S? (« = 1, 2). (19) 
The SO (10) symmetry breaking proceeds as follows: 

SO(10) x U(1) PQ ^ SU(3) C x SU(2) L x SU(2) fl x U(l) B _ L x U(1) PQ 

4 ^ 6 SU(3) C x SU(2) L x U(l) y 3 SU(3) C x U(l) em , (20) 

where the Higgs fields necessary to implement this chain are as indicated. Under U(1)pq, the 
Higgs fields transform as follows: 

(21O) ^ (21O) ; (126) ^ e 2i^(126) ) (45) ^ ^(45^ (1O) ^ e -«yi0)_ (21) 

As in the fermion case, these U(1)pq transformation properties ensure that the action of the 
residual PQ symmetry on these fields is identical to that of the center of SO (10). Note that 
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Figure 4: Gauge coupling unification in the SO(10) axion model at a scale Mx = 4.13 x 10 
GeV, with the addition of extra Higgs components at an intermediate scale Mb-l — 8-54 x 10 
GeV. Here a{ 1 = \oi7 B _ L y 2 + |«2 1 at scale Mb-l- 



all Higgs fields except for </>( 210 ) are complex. The allowed Yukawa couplings are (in schematic 
form) 

^160 (1O) , ?/Wl60 (126)t , ?/WlO0 {45) . (22) 

The allowed Higgs couplings include 

^(210)^(126)^(126)^(45)^ ^(210) 0(126)^(10) 0(45) ? 0(210)^,(126) 0(10) (33) 

These couplings guarantee that U(1)pq is the only global symmetry present. They also guaran- 
tee that 0*- 45 - 1 — > — 0^ 45 - ) is not a summetry of the Langragian so that the domain-wall problems 
associated with this symmetry can be avoided [27]. Now, (0^ 210 ^) cannot break U(1)pq since it 
is neutral under U(1)pq. Hence U(1)pq is broken at the intermediate scale Mb-l by ((// 126 )), 
(0( 45 )>. 

To make contact with the discussion on CW-based SU(5) inflation, we need to make sure 
that the tree-level couplings in Eqs. (j22j) and (|2^|) respect scale invariance. In particular, any 
cubic scalar coupling which respects SO(10) xU(l) pq is accompanied by the inflaton field 4> 
and suitable dimensionless coefficients. 

In its minimal form, the SO(10) axion model includes the following Higgs contributions in 
the renormalization group equations (RGEs): the (1,1,3,-1-1) and (1,3,1,-1) components of the 
126-plet, the (1,3,1,0) and (1,1,3,0) components of the 45-plet and the (1,2,2,0) components 
of the 10-plet are included between Mx and Mb-l, while only the SM Higgs doublet in the 
(1,2,2,0) components of the 10-plet is included between Mb-l and Mz- The fermions in the 
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10-plet acquire masses ~ M B _ L . Using two loop RGEs for the standard model gauge couplings 
I2HJ, we find M B _ L = 3.5 x 10 11 GeV and M x = 1.05 x 10 15 GeV. 

To rectify the situation and obtain energy scales ~ 10 16 GeV, we consider threshold correc- 
tions from suitable (intermediate mass ~ Mb-l) scalar components in the 45 and 126 multi- 
plets. To illustrate this, we consider two appropriate combinations of these scalar components. 
In order to simplify matters, we include only the left-right symmetric scalar components in each 
of these cases. For our first example, we take (8, 1, 1, 0) from 45, (3, 1, 1, a/3/2(— 1/6)) + (3, 1, 
1, y/s/2(l/6)) and (3, 1, 3, y/3/2(-l/6) + (3, 3, 1, y/s]2(l/6)) from 126 all at scale M B -l- Em- 
ploying the two loop RGEs for the gauge couplings, we obtain unification at M x = 2.61 x 10 16 
GeV with M B -l = 8.0 x 10 11 GeV. 

For our second example, we choose at scale Mb-l the following multiplets: (8, 1, 1, 0) and 
(3, 1, 1, y/3/2(-2/6)) + (3, 1, 1, y/3/2(2/6)) from 45, and (1, 2, 2, 0), (3, 1, 3, ^/3/2{-l/6) + 
(3, 3, 1, v /372(l/6)) and (3, 2, 2, y/t/2(-2/6) + (3, 2, 2, yjt/2{2/6)) from 126. In this case, 
we obtain unification at M x = 4.13 x 10 16 GeV with M B -l = 8.54 x 10 11 GeV (see Fig [p. 



Reheat Temperature and Non-Thermal Leptogenesis 

In SU(5), it is natural to consider the Yukawa coupling of the SU(5) singlet inflaton field to 
right-handed neutrinos Ni (i = 1, 2, 3). For simplicity we consider a simple Yukawa coupling 

Y N (f>NN, (24) 

where Y/v ~ / ((f)) is the coupling strength. For direct decay of the inflaton via this coupling, 
we take 2M/v < m^; then, using (0) ~ (10 - 100) mp, we expect Y/v ^ 10~ 6 . Thus, we estimate 
the maximum value of the reheat temperature T r to be 

T r ~ (^) K ^p) 1 ' 2 ~ 10 9 GeV. (25) 
In the SO(10) model, the value of T r is estimated from the coupling [22] 

126 x 126 x 2 , (26) 

which carries a negative sign in front in order to induce the intermediate scale breaking. (There 
is a similar coupling involving the 45-plet.) The coefficient of the expression in Eq. ( 1261) is of 
order (Mb-lI (0)) 2 - Similarly as above, we take 2Mb-l ^$ ^<j>- The decay rate of the inflaton 
into the scalar components (10,1,3) of the 126-plet is then given by 



Mb-lVm 2 



M J 

The reheat temperature T r is estimated to be 



(27) 



T r ~Vr^~(^)(Jf)AW (28) 
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The scalar bosons produced by the inflaton rapidly decay through their Yukawa couplings into 
right-handed neutrinos. Under the assumption of hierarchical right-handed neutrino masses, 
the lepton asymmetry is given by [23J [2U ED] 

ni/s < 3xl0 -.»(^)(^)(^), (29 , 

where M^. denotes the mass of the heaviest right-handed neutrino. From the experimental 
value of the baryon to photon ratio tjb = 6.225 x 10~ 10 [3TJ, the required lepton asymmetry is 
found to be n^/s ~ 2.5 x 10 -10 . With T r ~ (10 6 -10 8 ) GeV, these heavy neutrinos, with masses 
on the order of 8 x (10 9 -10 n ) GeV or so, can give rise to the observed baryon asymmetry via 
non-thermal leptogenesis. 

Magnetic Monopoles, Axion Domain Walls and Inflation 

The spontaneous breaking of any GUT symmetry yields topologically stable magnetic monopoles 
In our SU(5) and SO(10) examples, these monopoles carry masses of order where «g ~ ~ 
2jj denotes the GUT coupling constant. In the SU(5) model, for instance, monopoles are pro- 
duced during inflation once the magnitude of the coefficient ^((ft) 2 associated with the TrS 2 

term in Eq. (pQ) exceeds ~ H 2 , where H = ( ~ (°- 01 ~ 4 ) x 1Ql4 GeV ) denotes the 

Hubble constant during CW inflation. Using Tableland Eq. ([3]), it is readily checked that 
this occurs for values of 0o- Thus, the corresponding number of e- foldings is much greater 
than 50-60 and the monopoles are inflated away. This holds for the inflaton field rolling to its 
minimum from values smaller than its present VEV M. Clearly, the same results hold if the 
inflaton rolls to its minimum from values that are larger than M. 

The breaking of U(1)pq symmetry occurs at intermediate scales ~ 10 12 GeV <C H. In this 
case the corresponding phase transition is completed at the end of inflation, with T r ~ 10 6 -10 9 
GeV. Thus, it is prudent to ensure that there are no stable domain walls in the model by 
including in the SO(10) model two 10-plets of new matter fields (see Eq. (ITU1) ). 

GUT Inflation and Proton Decay 

Our analysis above based on GUT inflation shows that the vacuum energy during inflation, 
for both the CW and the Higgs potential, is of order 10 16 GeV if we wish to obtain a scalar 
spectral index close to 0.96, as suggested by WMAP5. In the CW case, the SU(5) and SO(10) 
superheavy gauge bosons which mediate proton decay, are then expected to possess a mass Mx 
which is a factor 2-4 larger than this scale. For the Higgs potential, the connection between 
the estimated scale of vacuum energy and M x is somewhat tenuous, but it seems reasonable 
to assert that they are of comparable magnitude. Put differently, the value for Mx inferred 
above is a factor ~ 1-10 larger than the lower bound ~ 4 x 10 15 GeV obtained from the 
observed stability of the proton [8]. Thus, we expect proton lifetime estimates based on gauge 
boson mediated decays to be more in line with what one finds in supersymmetric GUTs, where 
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Coleman- Weinberg Potential 


Higgs Potential 


M x ~ 2 Vo 1/4 (GeV) 


r(p — > 7r°e + ) (years) 


M x ~ U 1/4 (GeV) 


r(p — > 7r°e + ) (years) 


5.0 x 10 15 


1.8 x 10 34 


1.0 x 10 16 


2.8 x 10 35 


1.0 x 10"* 


2.8 x 10 35 


1.2 x 10 ie 


5.8 x 10 35 


1.2 x 10 16 


5.8 x 10 35 


1.4 x 10 16 


1.1 x 10 36 


1.8 x 10 16 


2.9 x 10 36 


1.6 x 10 16 


1.8 x 10 36 


2.2 x 10 16 


6.6 x 10 36 


1.8 x 10 16 


2.9 x 10 36 


2.7 x 10 16 


1.5 x 10 37 


2.1 x 10 16 


5.5 x 10 36 


3.5 x 10 16 


4.2 x 10 37 


2.4 x 10 16 


9.3 x 10 36 


6.0 x 10 16 


3.6 x 10 38 


2.9 x 10 16 


2.0 x 10 37 



Table 3: Superheavy gauge bosons masses and corresponding proton lifetimes with a G = ^ in 
the CW and Higgs models. Note that since the lifetime depends only on Mx, the results shown 
here apply equally well to the BV and AV branches in each model. 



Mx ~ 2 x 10 16 GeV [32]. A few lifetime estimates based on the naive expression for the decay 
rate [8] 

IV-^jjjjh (30) 

are presented in Table [3j (Here m p is the proton mass.) 

The presence of an intermediate scale ~ 10 12 GeV in the SO(10) axion model leads to the 
appearance of intermediate mass scalar leptoquarks (3, 1, —2/3) which mediate proton decay. 
In this case, the proton decay rate is given by [8] 

r P ^l^l 2 v74^, (31) 

B—L 

where Y u ^ denote the Yukawa couplings of the u and d quarks respectively. For an intermediate 
scale ~ (0.5-1) x 10 12 GeV as in Fig. HI the proton lifetime is estimated to be of order 10 34 -10 35 
years. 



Summary 

Assuming that inflation is associated with a GUT symmetry breaking phase transition (so that 
primordial GUT monopoles do not pose a cosmological problem) , we have provided estimates for 
the energy scale of observable inflation which yield values for the scalar spectral index centered 
close to 0.96. For both CW and Higgs inflationary potentials, an estimate for the masses of the 
superheavy gauge bosons which mediate proton decay can then be derived. Masses of order 
10 16 GeV for these bosons are favored, yielding proton lifetimes of order 10 34 -10 38 years. In 
some SO(10) models with axion dark matter, lifetimes of order 10 34 -10 35 years are realized 
with proton decay being mediated by scalar leptoquarks possessing intermediate scale masses 
of order 10 12 GeV. 
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